Generic strong coupling behavior of Cooper pairs in the surface of superfiuid nuclei 
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With realistic HFB calculations, using the D1S Gogny force, we reveal a generic behavior of 
concentration of small sized Cooper pairs (2-3 fm) in the surface of superfiuid nuclei. This study 
confirms and extends previous results given in the literature that use more schematic approaches. 
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The opportunities offered by the new radioactive beam 
facilities to study the properties of weakly bound nuclei 
with large neutron skins or halos triggered new interest 
for the issue of space correlations induced by the forma- 
tion of Cooper pairs. The spatial correlations of Cooper 
pairs in superfiuid nuclei have not been extensively stud- 
ied in the past, but nevertheless a certain number of in- 
vestigations, some rather early, do exist. Mostly this was 
done for the single Cooper pair problem. For example the 
rms diameter of the extra neutron pair in 18 O is shown as 
a function of the nuclear radius by Ibarra et al. [lj . One 
sees a strong minimum in the nuclear surface, indicating 
an rms separation between the two active neutrons of the 
order of 2-3 fm. A similar behaviour was found later by 
Catara et al. Q and Ferreira et al. [H for a neutron pair 
in 206 Pb and 210 Pb. More recently, there are also many 
investigations of the single Cooper pair problem in the 
halo nucleus 11 Li 

One of the rare papers where spatial correlations of 
Cooper pairs are investigated in superfiuid nuclei is the 
one of Tischler et al. @ where the probability distribu- 
tion of the pairs is shown as a function of the center of 
mass R— | |ri +r 2 1 and the relative distance of the nucle- 
ons in the pairs r — |r\ — r 2 | with (r\, f 2 ) the coordinates 
of two nucleons. They showed that in the open shell iso- 
tope 114 Sn one also finds Cooper pairs with short range 
space correlations, like in one pair systems. They con- 
firmed also the finding of Catara et al, i.e., the important 
role played by the parity mixing for inducing short range 
space correlations. Most of those older works were, how- 
ever, done using rather schematic models and/or pairing 
forces. There exists, however, one study with a realistic 
pairing force (i.e., the Gogny interaction) by Barranco et 
al. 0], dedicated to nuclei embedded in a neutron gas, a 
system found in the inner crust of neutron stars. One of 
the first systematic analyses of strong di-neutron spatial 
correlations induced by the pairing interaction was done 
recently by Matsuo et al.Q, using a zero range pairing 
force. The study of nuclear surface pairing properties 
was also the aim of several half infinite matter investiga- 
tions 0] , [l(| • It was found that the pair density reaches 
out further than the ordinary density but neither the lo- 
cal coherence length nor the probability distribution of 
the pairs were calculated. 

The aim of the present work is to verify how much all 



these relative scattered pieces of information withstand 
a general study of superfiuid nuclei using one of the 
most performant HFB approaches, that is employing the 
finite range Gogny DIS-interaction |12j| . As a matter 
of fact we will see that many of the earlier findings 
are qualitatively or even quantitatively confirmed. 
Indeed, we will show that the strong concentration of 
pair probability of small Cooper pairs in the nuclear 
surface is a quite general and generic feature and that 
nuclear pairing is much closer to the strong coupling 
regime [g, [l3| than previously assumed. 
We will start by explaining shortly how the spatial 
properties of nuclear pairing are investigated within the 
HFB approach and then we shall present our results 
and conclusions. For further understanding of the 
phenomena, a simple semiclassical analytic model for 
nuclear pairing will be also considered. 

It is well known [llj that pairing correlations can be 
adequately studied with the Cooper-pair probability \k\ 2 , 
where (in standard notation): 



,{n,r 2 ) =< HFB\^,{r 2 )Mn)\HFB > 



(1) 



is the anomalous density matrix or pairing tensor in r- 
space, calculated with the HFB-wave function. Following 
Refs. 0, [fjj], we shall also consider the quantity: 



P(R,r) = R 2 r 2 \n(R,rf 



(2) 



which is the pair probability averaged over the angle be- 
tween R and r and multiplied by the phase space factors 
R 2 r 2 . This quantity is important since it determines the 
two-particle spectroscopic factor |14| and other expecta- 
tion values of two-body operators (e.g., pairing energy). 
Let us mention that Eq.([2]) is formally the same for the 
single Cooper pair problem [l|, [1J!, Q and for the case 
of Cooper pairs in a condensate Q . 

In what follows, we shall consider the HFB expression of 
the pairing tensor k in center of mass and relative coor- 
dinates given by: 
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FIG. 1: 
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P(R, r) calculated with HFB-D1S for 22 O, 60 Ca, 60 Ni, 
120 Sn, 212 Pb. Scales have been multiplied by a factor 



where < nlNl; 0|mZiri2Zi; > is the Brody-Moshinski 
bracket, u n i(r) are the radial wave functions of the har- 
monic oscillator and kS u is the matrix of the pairing ten- 
sor for a given angular momentum Ij. As defined here, 
the latter has an intrinsic parity (— The HFB calcula- 
tions are performed with the D1S Gogny force [15J. In 
the calculations a basis with 15 harmonic oscillator shells 
have been considered. The contour-lines of P(R, r) for 
various superfluid nuclei are shown in FIGQ] The strik- 
ing feature is that for all these nuclei the same scenario, 
with only slight modulations, emerges: the probability 
P(R, r) is strongly concentrated in the surface with a 
small diameter of the pairs of the order of 2 — 3/m. 
In FIG[TJ we show nuclei close to the neutron drip- line 
( 60 Ca) as well as nuclei closer to stability. Seemingly, 
there is no essential difference in the behavior of P{R, r) 
between very neutron-rich nuclei and stable ones. This 
fact explains why one finds in all the superfluid nuclei a 
high probability for two-neutron transfer reactions. 

To conclude from the strong concentration of P(R, r) 
in the surface that Cooper pairs are mostly sitting in the 
surface of the nucleus may be, however, a bit mislead- 
ing. In FIG[2j we show the pair probability \k(R, r)\ 2 ) 
without the factor R 2 r 2 which enters in the probabil- 
ity P(R,r). One can notice that for all tin isotopes 
the Cooper pairs have very small extension in r-direction 



{J r 4 \ K (R,r)\ 2 dr)^ 2 
(f r 2 \ K {R,r)\ 2 dry/ 2 



(4) 



FIG. 2: \k(R, r)| 2 calculated with HFB-D1S for 104 Sn, 120 Sn, 
128 Sn. Scale has been multiplied by a factor of 10 6 . 



throughout the nuclear radius. However, the distribu- 
tion in R is rather different in the three isotopes. The 
difference comes from the localisation properties of the 
single-particle shells which are closest to the chemical po- 
tential. Thus, the pronounced concentration of \k(R, r)\ 2 
laround R=5 fm in 128 Sn is due to the surface localisa- 
tion of the single-particle wave function l/in/2, which 
becomes much closer to the chemical potential in this 
isotope compared to lighter ones. One can also notice 
| that in 120 ' 12S Sn the pair probability has a sizeable value 
for small values of R, which comes mainly from the con- 
tribution of the state 2>Si/2 to pairing correlations. If 
we have had chosen the neutron deficient Pb-isotopes, in 
which there is no s-state in the major shell, one would 
rather see a depression of pair probability at the origin. 
Therefore, to say where the Cooper pairs are preferen- 
tially located in nuclei is a somewhat subtle question be- 
cause the answer depends rather strongly on the shell 
structure (see also pif). The shell structure dependence 
of \k(R, r)\ 2 is largely washed out by the phase factor 
r 2 R 2 when P(R,r) is calculated. We shall make further 
investigations on this issue in a future work. However, we 
want to point out again that in all expectations values, 
like e.g. the pairing energy and spectroscopic factors, it 
is P(R,r) which counts and not \k(R, r)\ 2 . 

The fact that the Cooper pairs with small size are con- 
centrated in the surface can be also seen from the depen- 
dence of the coherence length on the center of mass of 
the pairs. The coherence length is defined as: 



It is shown for various nuclei in FIG(31 One sees well 
defined and pronounced minima at £ ~ 2 — 3/m for R 
of the order of the surface radius. As we have already 
mentioned, a small coherence length in the case of a 
single Cooper pair has already been found previously 
for ls O in Ref.[l|l. It is also the case for the Cooper 
pair in 11 Li Our calculations did not allow to 

go much beyond the minima because of the employed 
harmonic oscillator basis what becomes inaccurate far 
outside the nuclear radius. However, the position of 
the minima is always clearly identified and seen to be 
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FIG. 3: Coherence length for 22 O, b0 Ca, 

120 Sn, 212 Pb. 
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FIG. 4: Contributions of definite parity to P(R,r) calculated 
with HFB-D1S for 120 Sn. 



similar in all cases. What is surprising is that the size of 
the Cooper pairs starts to decrease already well inside, 
around R = 2/m. Moreover, the decrease towards the 
surface is approximately linear. 

In order to demonstrate that the strong concentration 
of small Cooper pairs in the surface of the nuclei is not 
a trivial effect, we decompose K,(R,r) in a part K e (R,r) 
which contains only even parity wave functions and a part 
K (R,r) which contains only the odd parity ones, i.e., 
n(R,r) = K e (R,r)+K (R,r). In FIG (H we show what are 
the probability distributions for P e (R,r), P a (R,r) and 
P eo (R,r) in the case of 120 Sn. The quantity P eo (R,r) 
corresponds to the interference term 2K e n . From FIG. 4 
one can see that selecting only either even or odd parity 
states in k(R, r) has a strong delocalisation effect on the 
Cooper pairs: they are democratically distributed with 
respect to an interchange of R and r variables (one should 
notice that in Eq.(3) the symmetry between R and r in- 
volves a factor 2, which comes through the standard defi- 
nition of Brody-Moshinsky transformation) . So no small 
Cooper pairs in the nuclear surface are prefered at all 
in those cases. The concentration only shows up, when 
even and odd parity states are mixed. This is clearly re- 
vealed in looking at the interference term P eo (R,r). We 
see that it is negative for regions close to the r-axis and 
positive close to the R-axis. We checked that this sce- 
nario stays the same for all other superfluid nuclei con- 
sidered. This scenario was also nicely described in the 
papers by Catara et al. @ and Tischler et al. Mix- 
ing of parities naturally occurs in heavy nuclei because of 
the presence of intruder states of unatural parity in the 
main shells of given parity. However, as seen in FIG[T] 
first panel, the concentration of Cooper pairs also occurs 
in such light nuclei as the oxygen isotopes where no in- 
truders are present. This means that pairing in nuclei 
is sufficiently strong so that k grabs contributions from 
several main shells, allowing for parity mixing even in 
light nuclei. If one artificially restricted the pairing con- 
figurations, e.g. in 22 0, to the s-d shell, then certainly no 
concentration effect at all would be seen (in this respect, 



see also the study in @). So to grasp the full physics 
of nuclear pairing it is very important to work in a large 
configuration space, comprising several shells below and 
above the active one (see also [8|). 
In order to understand in more detail where this extraor- 
dinary concentration effect from even-odd parity mix- 
ing comes from, let us consider a very simple model. 
We got inspired by the Thomas-Fermi model presented 
in Ref. where the anomalous density matrix is given 
by K TF (R,r) ~ kF{R)jo{kF{R)r) with jo(x) a spheri- 
cal Bcssel fonction, k F (R) = ^J^{fJ.-V(R)) 6(/x - V) 
the local Fermi momentum, \x the chemical potential (or 
Fermi energy), and V(R) is phenomenological mean field 
potential. It can be shown [18[ that a slightly more elab- 
orate semiclassical version can be written as 

771 f 

KSC ( R > r ) = Jp~2 J dE K ( E ) k E(R) 3o(k E (R)r) , (5) 

where fc^(i?) is the local momentum at energy E, ob- 
tainable from k F in replacing u by E, and k(E) is the 
continuum version of the k's for the individual quantum 
levels: k{E) = A(E)/(2 V / {E - a) 2 + A(E) 2 ). We see 
that for very small A's, one gets back the TF model 
However, for realistic gap values the distribution of k's 
is very important, otherwise the concentration effect will 
not show up. For A(E) we adjust a Fermi function to 
represent on average the gap- values of the individual sin- 
gle particle levels. An example can be seen on FIG[4] of 
Ref. [2(| • In the present work, we have fitted the function 
A(E) on HFB-D1S results for 120 Sn. A good fit function 
is given by A(E) = 4/[l + exp(E — /x)/20] (all numbers 
are in MeV). For the mean field potential V(R) we take 
the Woods-Saxon form of Ref. [13] • The chemical poten- 
tial jj, is determined, as usual, via the particle number 
condition. 

In FIG [5] we show the corresponding semiclassical prob- 
ability P sc (R,r). We see qualitatively good agreement 
with the quantal HFB-results, for instance in what con- 
cerns the concentration of small Cooper pairs in the sur- 
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FIG. 5: Parities contributions P sc (R,r) calculated with the 
semiclassical model for 120 Sn. 



face. We also show in FIG|5]the parity projected prob- 
abilities. As in FIGJH one sees the strong delocalisa- 
tion effect. In our model this can be understood analyt- 
ically. Parity projection can be written as K e / (r*i, r*2) = 

\Wi&) ± «(ri,-r 2 )] - f) ± «(f,2fl)] (see 

Ref. [131). We therefore see that good parity implies, 
up to a scale factor, a symmetrisation in coordinates R 
and r. This is general and can be investigated analyti- 
cally in the TF model. We also calculate the coherence 
length semiclassicaly. We find qualitatively the same be- 
havior as in the quantal calculation of FIGlH 
The analytic model also allows to quickly grasp the sig- 
nificance of the coordinates used by Matsuo et al. [8|. 
There, one takes a reference particle at position r\ on 
the z-axis, i.e. fi = zi& z . Moving the second particle on 
the z-axis we see that for P e / a two symmetric peaks at 
r*2 = z\e z and ?i — —z\e z appear whereas for the total 
probability only one peak on the side of the test parti- 
cle appears. This is a clear signature of strong pairing 
correlations as also pointed out in @. 



In conclusion, we showed that Cooper pairs in su- 
perfluid nuclei preferentially are located with small size 
(2 — 3/m) in the surface region. There, they maximally 
profit from the Cooper phenomenon, that is, with re- 
spect to the neutron-neutron virtual S-state in the vac- 
uum (rms 12 fm, 0]), strong extra binding occurs, as 
long as the density is not too high. Further to the center 
of the nucleus the stronger effect of the orthogonalisation 
of the pair with respect to the denser core-neutrons per- 
turbs the pair wave function. It starts to oscillate and 
expands again ||. That this simple, physically appeal- 
ing and generic picture, is so pronounced, has come as a 
surprise. It is certainly important for the interpretation 
of pair transfer reactions. Most of these facts had al- 
ready been revealed in the past for specific examples and 
schematic models and forces. We think, it is the merit of 
this work that it demonstrates with realistic HFB calcu- 
lations using the finite range D1S force the generic aspect 
of strong coupling features of singlet isovector pairing in 
nuclei. These features are in agreement with the ones 
recently put forward by Matsuo et al. Q. Let us men- 
tion that the strong coupling features revealed here are 
somewhat contrary to the old believe [22j that the coher- 
ence length of nuclear pairs is of the same order or larger 
than the nuclear diameter. On the contrary, a much more 
diverse local picture has emerged. This may also be the 
reason for the rather good succes of LDA for nuclear pair- 
ing found in the past [21| . In spite of the strong coupling 
aspects revealed in this work, we hesitate to say that 
there is Bose-Einstein condensation (BEC) of isovector 
Cooper pairs, since this, strictly speaking, occurs only 
for ( in infinite matter) negative chemical potential, what 
means true binding. However, fj, never turns negative for 
isovector pairing in infinite nuclear or neutron matters. 
Nuclear isovector pairing is just in the transition region 
from BEC to BCS. 

Acknowledgements : We would like to thank Marc 
Dupuis for his help in the HFB code. 



[1] R.H. Ibarra, N. Austern, M. Vallieres, D.H. Feng, 

Nucl.Phys. A288, (1977) 397. 
[2] F.Catara, A. Insolia, E. Maglione, A. Vitturi, Phys.Rev. 

C29 (1984) 1091. 
[3] L. Ferreira, R. Liotta, C.H. Dasso, R.A. Broglia, A. 

Winther, Nucl.Phys. A426 , (1984) 276. 
[4] G.F. Bertsch, H. Esbensen, Ann. Phys. 209 (1991) 327; 

H. Esbensen, G.F. Bertsch, K. Hencken, Phys. Rev. C56 

(1999) 3054. 

[5] K. Hagino, H. Sagawa, J. Carbonell, P. Schuck, submit- 
ted to publication arXiv:nucl-th/0611064vl. 

[6] MA. Tischler, A. Tonina, G.G. Dussel, Phys. Rev. C58 
(1998) 2591. 

[7] F. Barranco, R.A. Broglia, H. Esbensen, E.Vigezzi, Phys. 

ReV. C58 1998 1257. 
[8] M. Matsuo, K. Mizuyama, Y. Serizawa, Phys. Rev. C71 

(2005) 064326. 



[9] M. Farine, P. Schuck, Phys. Lett. B459 (1999) 444. 
[10] M. Baldo, U. Lombardo, E. Saperstein, M. Zverev, Phys. 

Lett. B459 (1999) 437. 
[11] P. Ring and P. Schuck, The nuclear many-body problem, 

Springer- Verlag (1980) 
[12] J.-F. Berger, M. Girod, D. Gogny, Comp. Phys. Comm. 

63 (1991) 365. 

[13] M. Baldo, U. Lombardo, P. Schuck, Phys. Rev. C52 
(1995) 975. 

[14] F. Catara, A. Insolia, E. Maglione, A. Vitturi, Phys. Lett. 
B149 (1984)41 

[15] J. Decharge, D. Gogny, Phys.Rev. C21, (1980) 1568. 
[16] N. Sandulescu, P. Schuck, X. Vinas, Phys. Rev. C71 

(2005) 054303. 
[17] S. Shlomo, Nucl.Phys. A539 (1992) 17. 
[18] X. Vinas et al., in preparation. 

[19] X. Vinas, P. Schuck, M. Farine, M. Centelles, Phys. Rev. 



■5 



C67 (2003) 054307. [22] A. Bohr, B.R. Mottelson, Nuclear Structure, Ben- 

[20] P. Schuck, K. Taruishi, Phys. Lett. B385 (1996) 12. jamin(1975). 
[21] H.Kucharek, P.Ring, P.Schuck, R.Bengtsson, M.Girod, 

Phys.Lett. B216 (1989) 249. 



